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ABSTRACT 

Molecular line observations may serve as diagnostics of the degree to which the num- 
ber density of cosmic ray protons, having energies of 10s to 100s of MeVs each, is 
enhanced in starburst galaxies and galaxies with active nuclei. Results, obtained with 
the UCL_PDR code, for the fractional abundances of molecules as functions of the 
cosmic-ray induced ionisation rate, £, are presented. The aim is not to model any par- 
ticular external galaxies. Rather, it is to identify characteristics of the dependencies of 
molecular abundances on £, in part to enable the development of suitable observational 
programmes for cosmic ray dominated regions (CRDRs) which will then stimulate de- 
tailed modelling. For a number density of hydrogen nuclei of of 10 4 cm -3 , and high 
visual extinction, the fractional abundances of some species increase as £ increases to 
10~ 14 s _1 , but for much higher values of C the fractional abundances of all molecu- 
lar species are significantly below their peak values. We show in particular that OH, 
H2O, Ht , H30 + and OH + attain large fractional abundances 10 -8 ) for ( as large 
as 10 -1 ^ s -1 . HCO+ is a poor tracer of CRDRs when ( > 10 " 
species may be useful tracers of CRDRs gas in which £ ~ 10 
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Sulphur-bearing 
. Ammonia has 

a large fractional abundance for £ ^ 10~ ie s _1 and nitrogen appears in CN-bearing 
species at significant levels as £ increases, even up to ~ 10~ 14 s -1 . In this paper, we 
also discuss our model predictions, comparing them to recent detections in both galac- 
tic and extragalactic sources. We show that they agree well, to a first approximation, 
with the observational constraints. 
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INTRODUCTION 

The high spatial density of massive star formation in mergers 
and s t arburst galaxies (e . g . lAcero et al.ll2009l ; ISuchkov et al.l 
ll993l ; lAcciari et al.ll2009h creates regions of extremely high 
cosmic ray energy density, up to about ten thousand t imes 
that in the Milky Way Galaxy. iPapadopoulosI (2010) has 
proposed that these large energy densities alter the heating 
rates and ionization fractions in dense gas (n(H~2) > 1 x 10 4 
cm' 3 ) in the UV-shielded cores that contain much of the 
molecular gas in these galaxies, so that these cosmic ray- 
dominated regions (CRDRs) have different initial conditions 
for st ar formation l|Papadopoulos|[201ol : iPapadopoulos et"al] 
2010). These conditions affect the subsequent evolution of 
the gas and may even lead to a t op-heavy Initial Mass Func - 
tion and bimodal star formation |Papadopoulos et alfc oiO). 
It is therefore important to determine whether there 



may b e useful molecular tracers of CRDRs. IPapadopoulosI 
I2OIOI ) has suggested several possible chemical signatures, 
but it seems necessary to make a fairly complete and self- 
consistent thermal/chemical model of dense gas subjected 
to very high fluxes of cosmic rays. Historically, the cosmic 
ray ionization rate in Galactic diffuse and molecular clouds 
has normally been determined by treating it as a free pa- 
rameter an d matching the predicted chemistry to that ob - 



serv ed (e.g. iBlack fc Dalgarnol Il977l : lHartquist et all Il97g| a 



* E-mail: eb@star.ucl.ac.uk; 
twh@ast.lceds.ac.uk; sv@star.ucl.ac.uk 



daw@star.ucl.ac.uk; 



and l Hartauketal l l ia lSnierrTt is necessary to r=e 
that procedure and to compute the va riation of the chem- 
istry as the ionization rate is increased. iLepp fc Tinl (|l998l ) 
made such a calculation and computed the chemistry in 
dense gas subjected to ionization rates varying over a wide 
range. Ho wever, their calculation was for a fixed tempera- 
ture. L ater, IMeiierink fc Spaansl (120 05). Spa ans k, Meiierinkl 
(2005), and IMeiierink et all (|2006t ) concentrated on a re- 
stricted chemical network with enhanced ionization rate. 
Recently. iBavet et al.l (|2009l ) compared the computed chem- 
istry for a large network of 131 species connected in 1700 
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reactions at just two particular cosmic ray ionization rates 
(the canonical Milky Way value and 100 times this value) in 
a self-consistent thermal/chemical and time-dependent one- 
dimensional mo del. In the present work, we use the same 
model as that of lBavet et alj (|2009l ) and compute the chem- 
istry for many ionization rates from Milky Way values up to 
about one million times larger. 

The thermal/chemical model is described in Section [2l 
and the input parameters are discussed in Section The 
outputs of the computations are the temperature and the 
chemical abundances as functions of visual extinction (from 
A v = 3 to 20 magnitudes) and of the cosmic ray ionisation 
rate (from ( = 2x 10~ 17 s _1 to ( = 5 x 10" 11 s _1 ). The 
model outputs are presented in Section [4] The relevance of 
our modelling to observations of galaxies is considered in 
Section [SI where we also give our conclusions. 



of species on to grain surfaces is excluded. The UCL_PDR 
code has be en validated again st all other commonly used 
PDR codes ijRollig et alJl2007t ). 

The model should be appropriate for regions with tem- 
peratures up to several thousand degrees. However, it is 
not valid for regions with temperatures approaching 10 4 
K and the associated high fractional ionisations. Thus, we 
present results only for models for which C is 10~ 12 s" 1 or 
less. Codes such as MOCASS IN ijErcolano et all 120051 ) or 
CLOUDY ijFerland et al.lll998h might be more suitable for 
models with higher values of The model does not include 
doubly ionised species other than Ca ++ , though cosmic- 
ray induced ionis ation will produce a variety of them (see 
IShaw et al.l 120081 ) . Charge transfer with neutral hydrogen 
atoms will contribute substantially to the removal of many 
of them in regions with £ ranging from 10 -14 to 10~ 12 s _1 . 



2 MODEL DESCRIPTION 



3 PARAMETER SELECTION 



The model used for our study of the che mistry in CRDRs 
is th e UCL_PDR model as described in iBell et al.l l|2005l . 
2006). This is a time-dependent Photon-Dominated Region 
(PDR) model with constant density. In the present appli- 
cation, the code is run for 10 s yrs by which time chemical 
steady state is reached in all models. Checks that we have 
made confirm that the fractional abundances at 10 6 yrs do 
not differ significantly from those at 10 8 yrs. However, even 
for high values of the chemistry requires close to 10 6 yrs 
to approach steady state. Though the timescale to convert 
atomic carbon to CO decreases with increasing (, neutral- 
neutral reactions limit the rates at which many of the ob- 
servable species form. Those rates do not depend sensitively 



on for values in excess of about 10 



However, they 



do decrease with decreasing metallicity. 

The code operates in one space dimension for an 
assumed semi-infinite slab geometry, and computes self- 
consistently the chemistry and the temperature as functions 
of depth and time within a semi-infinite slab, taking ac- 
count of a wide range of heating and cooling processes. In 
the present work, the code is used to determine the chem- 
ical and thermal properties at all depths up to 20 visual 
magnitudes. For t he present work, and consistent with the 
iBavet et al l (|2010l ) version of the code, additional radiative 
cooling due to rotational transitions of 13 CO, C 18 0, CS, OH 
and H2O h ave been implemente d. The escape probability 
formalism of Ide Jong et alj (|l980l ) is used to determine non- 
LTE level populations and resulting line intensities at each 
depth- and time-step, in the same manner as for the exist- 
ing coolants in the code. Collisional rates with H2 are taken 
from the Leiden At omic and Molecular Database (LAMDA; 
ISchoier et al]|2005h . The inclusion of these coolants allows 
the thermal balance to be more accurately determined in 
warm dense gas at high extinction. The UCL_PDR code 
also includes cooling due to the OI metastable levels, and 
the O - H + collisions and O - e~ collisions which excite 
them; these are important for the models for high values 
of The treatment of H2 self-shielding has also been up- 
da ted to use the r esults of detaile d calculation s perfo rmed 
bv lLee et al] l|l996h . as described in lBavet et all (|2010l ). The 
chemical network links 131 species in over 1834 gas-phase re- 
actions; only H2 is formed by surface chemistry; freeze-out 



The model requires the setting of a number of physical and 
chemical parameters. The choices we have made are listed 
in Table [1] We comment here on some of these. All the 
calculations are made for an assumed gas number density of 
1 x 10 4 hydrogen nuclei (in all forms) cm -3 . 

The cosmic ray ionisation rate is a free parameter in 
these calculations, and allowed to take values that range 
from C = 2 x 10~ 17 s" 1 to 1 x 10~ 12 s" 1 . 

The lowest value of £ selected here is comparable to 
that inferred for the Milky Way, while the highest value 
of £ is comparable to the highest values suggested to be 
appropria te for compact extr agalactic starburst regions (cf. 
eqn. 6 of lPapadopouioslbOlOt ). 

Galaxies are observed to have a r ange of metallicities 
fe.g. lBavet et al.ll2006l . IBell et al.ll2006t) ; we have considered 
a range from 0.1 to 4 times solar metallicity. Most of our 
computations are made using metallicity of 0.1 times solar, a 
value that may represent conditions in fairly young galaxies. 
However measurements of metallicity in some galaxies show 
both near-solar and larger than solar val ues (NGC 253, M 
82, IC 342 - see e.g. IZaritskv et alj 119941 ) . so we have also 
made calculations with higher metallicities. 

Several other parameters are assumed to scale linearly 
with metallicity. These are the dust:gas mass ratio, the H2 
formation rate, and the initial elemental abundances. The 
values of these parameters for solar metallicity are given 
in Table Q] The elemental abu n dance s given there are de - 
rived from Sembach & Savage il 19961) ; ISofia et al.l il 19971); 



(|l998l ); ISnow et al] (|2002l) : and iKnauth et alj 



I Meyer et al. 
((200j. 

The models require the specification of the mean dust 
grain radius and albedo, which affect the transfer of exter- 
nal radiation into the cloud. The values given in Table [1] are 
canonical. The external radiation field impinging on the one- 
dimensional semi-infinite slab is assumed to be 1000 Habing, 
much larger than for the Milky Way but possibly appro- 
priate for regions of active star formation in which CRDRs 
may be located. This value is comparable to values obtained 
fro m molecular line in ten sity modelli ng for several galaxies 
fe.g. lBavet et al.ll2006l and lAaltoll2008l ). The chemistry in the 
one-dimensional semi-infinite slab is explored as a function 
of depth or visual extinction, A v , from the edge (A v = 0) up 
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to A v = 20 mag. The model also requires for line-width cal- 
culations the specification of the microturbulen ce velocity. 
We use the value adopted bv lBavet etaD (|2009h . 



4 RESULTS 

We present in Fig [TJJ results for the temperature and 
chemistry at the three representative depths into the one- 
dimensional slab: A v = 3, 8 and 20 mag. 

The case A v = 3 is intended to represent translucent 
material, A v = 8 to represent the extended PDR gas (as 
detecte d in the nuclei of M 82, NGC 253, IC 342 and NGC 
4038 bv lBavet etafl 120091 ). while A v = 20 represents den se 
PDR components detected in galaxies (|Bavet et alj |2008l. 

Fig. [T] shows the results for temperature obtained from 
the self-consistent thermal balance, as a function of cosmic 
ray ionisation rate, at these three values of A v . All three 
curves show the same general behaviour, rising from low 
values at low ionisation rates to temperatures approaching 
1 x 10 4 K for the highest ionisation rates. As will be evi- 
dent in the plots of the chemical abundances, reactions at 
these very high temperatures destroy much of the chemistry. 
There is some difference in temperature between the three 
curves shown for low values of the ionisation rate. At A v = 3 
mag, the temperature lies between 50 and 100 K for ioni- 
sation rates less than about 1 x 10~ 14 s -1 . This relatively 
high temperature is maintained by the intense external ra- 
diation field. However at A v = 8 mag, the temperature falls 
to ~10 K at the lowest ionisation rates, as expected, since 
the external radiation field no longer plays a role. 

We note that the models give slightly higher temper- 
atures at A v = 20 mag than at A v = 8 mag. This is due 
to the fact that at high visual extinctions the C and CO 
lines become optically thick and hence are less able to cool 
the gas whilst the cosmic ray heating remains high and con- 
stant. For some models, there were some small difficulties 
associated with the convergence of the thermal balance. For 
example one sees small spikes in the thermal structure curve 
for A v = 3 and corresponding features in the plots of the 
fractio nal abundances. These dif ficulties are not due to bista- 
bility l|Boger fc Sternberg! I200r3 ). Rather they arise due to 
the stiffness of the coupled chemical and thermal balance 
equations. Their influence on the calculated abundances of 
observable species is not significant. The effects of high heat- 
ing rates on the chemistry is further investigated in Bayet 
et al. (2010). 

Fig-[U[3]and|3]present the results giving the abundances 
of various relevant species at the three specified values of A v . 
These species were selected from the 131 available partly be- 
cause many of them have already been detected in external 
galaxies, and partly to illustrate the sensitivity of species to 
the ionisation rate for a wide chemical variety. The range 
of abundances shown extends somewhat beyond the values 
that may be detectable. 

Although there are significant differences, the three 
cases are broadly similar in their behaviour with respect 
to the enhancement of the ionisation rate. As the ionisa- 
tion rate is increased from a value appropriate for the Milky 
Way, a rich chemistry is maintained up to a critical value of 
the ionisation rate, beyond which the molecular abundances 
fall rapidly as the ionisation rate is further increased. If the 



ionisation rate is increased to lxlO -12 s _1 , almost all the 
chemistry is effectively suppressed. The driver of this decline 
in molecular abundances is linked to the decline in molecu- 
lar hydrogen; this occurs when the ionisation rate is about 
lxl0~ 14 s" 1 , and H2 becomes a minor species when the 
ionisation rate is as large as 10 -12 s _1 . In these conditions, 
conventional astrochemistry - based on reactions with H2 
molecules - ceases. In this range of ionisation rate, the tem- 
perature rises abruptly to some thousands of Kelvin (see Fig. 
1), and hot atomic hydrogen is very destructive of molecules. 
However, we can see from the figures that potential molec- 
ular tracers can be identified for ionisation rates of about 
lxlO -13 s _1 but the molecular abundances drop substan- 
tially as £ increases much more. 

We now confine our remarks to the chemically richer 
regions of parameter space. There are some differences be- 
tween the panels in Fig. [3] (or [4}. Some species decline with 
increasing ionisation rate more rapidly than others. For ex- 
ample, sulphur- bearing species may be useful tracers of gas 
in which < = 10" 16 s" 1 but not when £ = 10" 14 s~\ Am- 
monia has a large fractional abundance for lower ionisation 
rates, but the abundance declines rapidly as the ionisation 
rate increases. Nitrogen appears in CN-bearing species at 
levels that are probably significant as the ionisation rate in- 
creases, even up to ~ 10 -14 s . 

There are some differences between the figures shown. 
At A v = 3, molecular abundances are significantly lower 
than those for higher visual extinctions. For example, CS is 
~ 10 4 times more abundant at A v = 8 than at A v = 3. The 
intense radiation field assumed in these models to impinge 
on the slab is responsible for inhibiting the chemistry at A v 
= 3. However, the model results for A v = 8 and 20 magni- 
tudes are essentially identical, as the external radiation field 
does not penetrate effectively to these depths. 

The behaviour of oxygen and carbon hydrides and their 
ions is somewhat different to the above behaviours. Some of 
these species show fractional abundances of rather high lev- 
els for ionisation rates of 10 -14 s , while others even sustain 
these high values for ionisation rates approaching ~ 10~ 12 
s _1 . For example, OH and H2O fractional abundances are 
~10~ r for an ionisation rate of ~ 10~ 14 s _1 , while OH and 
OH + may be as large as ~10 -8 even for an ionisation rate of 
~ 10 -12 s . The ions H^ and HsO" 1 " are likely to be abun- 
dant for £ = 10~ 13 s _1 ; however, the commonly used tracer 
HCO + is probably useless for values of £ larger than 10~ 13 
s _1 . The behaviour of these oxygen and carbon species is a 
consequence of the thermal stimulation of the endothermic 
reactions that may initiate oxygen and carbon chemistry by 
the increase in kinetic temperature arising from the higher 
ionisation rates. However, ultimately, the reduction in the 
H2 fraction suppresses these reaction networks. 

Given that the metallicity in some galaxies has been 
measured (see Section [5} and the values found to lie in a 
range often below but sometimes exceeding the solar value, 
it is worthwhile to explore the predictions of molecular abun- 
dances in our model for varying metallicity. Table [2] shows 
the computed fractional abundances for some atoms and 
molecules of observational interest for three values of the 
metallicity (0.1, 1.0 and 4.0 times the solar value), for both 
"high" (10 -13 s _1 ) and "low" (10 -16 s _1 ) cosmic ray ionisa- 
tion rates, and for two values of the visual extinction (A v = 
3 and 20 magnitudes). The results for the case when A v = 
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8 mag are closely similar to those for A v = 20 mag, and are 
not shown. The "high" io nisation rate is inten ded to repre- 
sent the case predicted bv lPapadopoulosI (|2010h for CRDRs, 
while the "low" case is expected to be roughly appropriate 
for galaxies similar to the Milky Way. 

As has been found by earlier studies of the dependence 
of chemistry on the metallicity, the behaviour is rather com- 
plicated. At first sight one may expect a higher metallicity to 
lead to larger fractional abundances of molecules. However, 
because of the chemical rate equations are strongly coupled, 
some species tend to increase with metallicity, while others 
decrease, for a given ionisation rate. Further, Table [2] shows 
that different behaviours may occur when the ionisation rate 
is changed. For example, the fractional abundance of C2H 
at A v = 3 mag with a high ionisation rate is a strongly in- 
creasing function of metallicity, while for A v = 20 mag and 
a low ionisation rate the C2H fractional abundance declines 
rapidly with increasing metallicity. 

CO is generally fairly strongly dependent on metallicity 
as it is formed in successive reactions. However, at A v = 3 
mag with a low metallicity, the molecule does not attain an 
abundance high enough to provide self-shielding against the 
radiation field which destroys this molecule through line ab- 
sorption. Thus, at low metallicity, the fractional abundance 
of CO is also low for A v = 3 mag. Consequently, HCO + is 
similarly affected. 

These kinds of behaviour can be understood in terms 
of the rates of the various competing processes. At A v = 
3 mag, the intense radiation field adopted in these calcu- 
lations (1000 Habing) may lead to photodissociation being 
a competitive process with ion-molecule reactions, while at 
higher values of A v photodissociation is negligible. The ion- 
isation rate drives the chemistry when photo-processes do 
not play a role, and high ionisation rates drive the chem- 
istry much faster than low rates. Thus, in the case of C2H, 
the dependence on metallicity is strong because it is formed 
by successive reactions involving C + , itself dependent on 
metallicity. Thus, at A v = 3 mag both photons and high 
fluxes of cosmic rays drive the chemistry strongly to create 
large abundances at high metallicities. However, at A v = 20 
mag with low ionisation rates, the C + and C2H abundances 
never attain high values, even at the highest metallicity. It 
is evident, therefore, that the choice of potential tracers de- 
pends sensitively on ionisation rate as well as metallicity. 



5 DISCUSSION AND CONCLUSIONS 

In this section we briefly explore the relevance of our theoret- 
ical predictions for galactic and extragalactic environments. 
We remind the reader that our attempt in this paper is not 
to model any particular source but rather to indicate qual- 
itatively if our model predictions (i.e. the variations of the 
species fractional abundances with respect to the changes 
in cosmic ray ionisation rate) agree to a first approximation 
with observations in the local and nearby Universe. 

Recent results indicate that higher cosmic ray ioniza- 
tion rates might be present in diffus e clouds throughout 
the Galaxy and i n the Orion bar (e.g. Ilndriolo et al.ll2007l ; 
IShaw et~al 1 120081 : Ilndriolo et all 120091 and even more re- 
cent results from H erschel/HIFI such as | Gupta et al.ll2010l ; 

(|2007h 



presents Hg~ column densities varying between 0.2— 2.1 x 10 1 
cm -2 depending on the diffuse cloud observed (see their Ta- 
bles 3 and 4). They derive estimates of £ ~ 0.5 — 3.2 x 
10 _16 s _1 . If we convert these column densities into frac- 
tional abundance estimates (i.e. 4.4 x 10" 9 — 4.4 x 10 -8 ) 
and compare them with those provided by our model show- 
ing the smallest A v (i.e. Fig. [2] for A v = 3 mag) and 
( = 2 x 10~ 16 s _1 , we find a very good agreement. Despite 
this encouraging agreement, the interpretation of the Hjj" 
observations is a highly debated subject and a more de- 
tailed analysis, requiring the modelling of a particular set of 
galactic sources (outside the sc ope of this paper) is r equired 
for any definitive conclusions. Ilndriolo et al.l (|2009h invoke 
a low-energy (~ 10 MeV) cosmic ray flux for a plausible 
explanation for the observed Hj results. Chemical models 
such as ours cannot be used to infer details of the cosmic 
ray spectrum responsible for inducing the ionization. 

Towards W49N. lNeufeld et all (|201Ch estimated column 
densities of OH + varying between 2.2 x 10 13 cm -2 and 
2.6 x 10 14 cm -2 whereas column densities of H20 + range 
from 4.2 x 10 12 cm" 2 to 2.7 x 10 13 cm -2 . These authors also 
inferred that these species are located in clouds of low molec- 
ular fraction. From their column densities we derive frac- 
tional abundances estimates of OH + of about 4.5 x 10~ 9 — 
5.4 x 10" 8 and for H 2 0+ of about 8.8 x 10" 10 - 5.6 x 10~ 9 . 
OH + fractional abundances are in agreement with our mod- 
els having £ > 5 x 10 _14 s _1 (see Fig. [2)| whereas H20 + 
fractional abundances are closer to predictions from models 
having £ ^ 10 — s — 1 (see also Fig. For values of £ in 
this range, Figure [2] also shows that n(H)/n(H2) becomes 
large, so that the ga s is mainly at o mic. This is consistent 
with the inference of lNeufeld et al.l d2010h t hat the gas has 
a low molecular fraction. Gerin et alJ ~ 201oT) detected OH + , 
H20 + and H3O 4 " towards the massive star-forming region 
G10.6-0.4 and estimated column densities of 2.5 x 10 14 cm -2 , 
6.0 x 10 13 cm" 2 and 4.0 x 10 13 cm" 2 , respectively. The OH+ 
and H2 O" 1 " column densities are roughly similar to those 
found bv lNeufeld et all (|2010h in W49N, which suggests that 
the HsO + fractional abundance in G10.6-0.4 may be of the 
order of ~ 10~ 9 . Figure [2] shows that values of this order 
can be achieved when £ is ~ 2 x 10" 14 s" 1 . 



Towards Orion-KL iGupta et all (|2010h showed that de- 
spite lower OH + and H20 + column densities than in the 
W49N case, a high cosmic ray ionisation rate (i.e. £ ^ 1 — 
2 x 10" 14 s" 1 ) is still required for reproducing these observa- 
tional values. Here, we found that models with £ ~ 10~ 14 s~ 1 
for OH + and models with ( 2x 10~ 14 s _1 agree with those 
observational column densities. These results als o seem con- 
sisten t, to a first approximation, with those of Sha w et al.l 
(2008) who showed that £ > 40 x the standard value (i.e. 



C ~ 7 x 10" 



is required to reproduce the abundances 



of several species in £ Persei. 

Interstellar molecular oxygen has been a target of 
SWAS, Odin, and - more recently - Herschel HIFI (See the 
talk given by P. Goldsmith at the Stormy Cosmos meeting: 



http://www.ipac.caltech.edu/ism2010/talks/Goldsmith 



iNeufeld et ail [2010). More precisely, Ilndriolo et al 



_StormyCosmos2010.pdf). The Herschel results confirm and 
extend the earlier results. Most galactic sources show no 
detectable emission, suggesting that in regions of modest 
temperature the O2 fractional abundance is low, with limits 
between a few times 10~ 9 to a few times 10~ 8 . Model 
results reported here show that at relatively low A v (Fig. 
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[2]) the O2 fractional abundance peaks at ~ 10 9 when 
C ~ 1 x 10~ 14 s _1 . At larger A v , then X(0 2 ) is ~ 10~ 8 when 
( ~ 1 x 10~~ 16 s _1 , but falls rapidly when £ is larger than 
~ 1 x 10 _14 s _1 . The model results appear to be consistent 
with the observational data. 

We consider then in particular for the extragalactic 
comparison, galaxies where AGN - and/or starburst - ac- 
tivities may have enhanced the cosmic ray ioni zation rates 
as we ll as, in some cases, created ULIRGs. In iBavet et al.l 

(2009) we selected well-known galaxies, such as Arp 220 or 
M 82 as examples of sources with active nuclei, and where 
therefore the cosmic ray ionization rate may be enhanced. 
Arp 220 is the prototypical ultraluminou s galaxy while M82 
is th e prototypical starburst. Recently (|van der Tak et ahl 
2008) H30 + has been discovered in these two regions and 
its fractional abundance has been estimated to be ~ 2- 
10 x 10 -9 . These authors found that observations of M 82 are 
matched by a high-£ PDR, i.e., an evolved starburst while 
X-ray models are best at reproducing the observations of 
Arp 220. In fact, our models indicate that one can obtain 
high abundance of this ion at low (i.e in a PDR) as well as 
high (i.e dense star forming gas) extinction as long as the C, is 
~ 10~ 13 s -1 and the metallicity is solar. At high extinction, 
which could represent the nuclear part of the galaxy, the 
abundance is higher, providing possibly a better match for 
the observations. We note here t hat we are not attempting 
to model Arp 220 that, since as Ivan der Tak et ail (|2008l ) 
pointed out, it has a quite unusual geometry. 

Another interesting object which has been recently 
studied in molecular emission is Mrk 231, a ULIRG. A 
high resolution SPIRE FTS spectrum r eveals the presence 
of ion s such as OH + , CH + and H20 + (|van der Werf et al.l 
2010). While abundances are not derived, we can use our 
Table [2] to determine what type of model is able to produce 
high fractional abundances 10 -10 ) of these three ions. 
We find that the only regime is an environment with low 
metallicity (0.1 solar), very high cosmic ray ionisation rates 
(> 10 ~ 16 s _1 ) and low visual extinction. Ivan der Werf et al.l 

(2010) explained the high abundances of these ions by in- 
volving XDR- c hemis try. We find that, in agreement with 
iPapadopoulosI (|2010l ). determining the origin of molecular 
emission from ULIRGs such as Mrk 231 is not trivial when 
both sources of energy (CR and X-rays) are present. 

Finally, our results are necessarily indicative rather than 
specific in that we do not estimate molecular line intensities. 
Our study was motivated by the recent investigation of fila- 
ments around the central galaxies of clusters of galaxies by 
IBavet et all l|2010l) and by the considerations of the effects of 
high c osmic ray ionization rates in ULIRGs bv lPapadopoulosl 
(2010) as well as recent Herschel results. We find that sev- 
eral species, many detected in extragalactic environments, 
are in fact tracers of very high ionization fractions. 

The general conclusions that we can draw from this 
study are as follows: 
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Figure 1. Temperature (in K) as a function of f (in s" 1 ) for 
A v = 3 mag (top), 8 mag (middle) and 20 mag (bottom). See the 
text in Sect. H] 



sation rates are enhanced, even up to rates as large as 10 1 
s" 1 ; 

• for galaxies with ionisation rates as large as 10 12 s 1 , 
gas of number density of 10 4 cm -3 is largely neutral and 
atomic, with a minor component of ions. Potential tracers 
of such gas are rare, but include molecular tracers OH and 
OH + and atomic tracers C and C + ; 

• Model results for the abundances of various molecular 
ions detected recently in the Milky Way and in external 
galaxies appear to be reasonably consistent with an origin 
in regions of enhanced cosmic ray flux; 

• if CRDRs are also regions of intense radiation fields, 
then translucent regions are chemically poor relative to re- 
gions that are more optically thick; 

• for regions with A v > 8 mags, the chemistry appears to 
be independent of depth; 

• the non-linear dependence of chemistry on metallicity 
may allow determinations of metallicity through molecular 
observations. 



• The dense gas kinetic temperature in galaxies with high 
cosmic ray ionisation rates is elevated above Milky Way val- 
ues. It reaches values ^ 3000 K for ionisation rates of about 
10- 12 s" 1 ; 

• many potential molecular tracers can be identified for 
dense gas in external galaxies in which the cosmic ray ioni- 



6 Bayet et al. 

Table 1. Model parameters (see text in Sect. 15) . 



Gas Density 




10 4 cm d 


Cosmic ray ionisation rate C, 


5xlCT 


11 T n 17 1 

11 -2xl0 Ll s 1 


Metallicity 




0.1, 1 and 4 zq 


Gas: dust mass ratio" 




100 


H2 formation rate a 


3xlO- 18 v / Texp( T5 |j)n 


(H)n H cm _3 s _1 


C/H a 


1.42 xl0~ 4 


0/H a 




3.20 Xl0~ 4 


S/H a 




1.43 xlO 


N/H a 




6.50 xl0~ 5 


Mg/H a 




5.10 xlO -6 


He/H a 




7.50 xlO" 2 


Si/H a 




8.21 xl0~ 7 


Fe/H a 




3.60 xl0~ 7 


Cl/H a 




1.10 xl0~ 7 


Na/H a 




8.84 xl0~ 7 


Ca/H a 




5.72 xlO" 10 


Grain size 




0.1 urn 


Grain albedo 




0.7 


External FUV radiation 




1000 Habing 6 


intensity 






A v maximum 




20 mag 


Microturbulence velocity 




1.5 krns" 1 



a : z = 1 = zq corresponds to solar values of the initial elemental 
abundance ratios while z = 1/10 zq means that the solar values 
of the initial elemental abundance ratios have been all divided 
by the same factor (of 10 in this example). The init i al elemental 
abund a nce ratios are f romlSembach fc Savage] (1996); Sofia et al. 



1997 1 ) : iMever et ail l|l99Sft : ISnow et al] l|2002l ); iKnauth et al 



2003). Similar assumptions are supposed for the H2 forma- 
tion rate coefficient and the dust:gas mass ratio; b : The unit of 
the standard Interstellar Radiation Field (ISRF) intensity is a 
mean inten sity of 1.6xl0 ~ 3 erg cm - 2 s _1 integrated over 912-2400 
angstroms jHabing|ll96sT) . 



Table 2. Fractional abundances (n(X)/nn) of species X (see col. 1), where nu is the total number of hydrogen atoms, for various metallicities and as a function of high cosmic ray 
ionisation rate (( = 1 X 10 — 13 s — x ) and lower cosmic ray ionisation rate (( = 1 X 10 — 16 s _1 ). The fractional abundances are given as a(b) which represents axlO b . The values listed are 
those obtained for opacities of A v = 3 mag and A v = 20 mag (see text in Sect. |4j. 



A v =3 mag 












A v =20 mag 












C (s" 1 ) 10" 13 


lO" 13 


lO" 13 


10 -16 


10 -16 


10 -16 


lO" 13 


lO" 13 


10~ 13 


10 -16 


10 -16 


10 -16 


metallicity (zq) 0.1 


1 


4 


0.1 


1 


4 


0.1 


1 


4 


0.1 


1 


4 


Molecule 



H 2 
OH 
H 2 
H 2 0+ 
C+ 

c 

CO 
CH+ 

CH 
CO+ 
C 2 H 
NH 3 

cs 
so 

H 2 CS 
H 2 S 
H s O+ 
HCO+ 
CN 
HCN 
HNC 
OH+ 
HOC+ 
H 2 CO 



2.69(-3) 
l.ll(-7) 
2.64(-10) 
2.03(-10) 
9.43(-6) 
4.76(-6) 
7.16(-9) 
3.82(-10) 
9.60(-10) 
8.12(-13) 
5.74(-16) 
9.31(-17) 
2.98(-17) 
1.68(-16) 
7.17(-24) 
2.00(-18) 
1.90(-13) 
1.73(-13) 
6.21(-12) 
6.85(-14) 
6.97(-18) 
2.52(-8) 
8.76(-14) 
3.53(-19) 



3.35(-l) 
1.19(-6) 
1.75(-7) 
4.25(-9) 
9.95(-5) 
3.06(-5) 
1.19(-5) 
6.37(-12) 
1.54(-9) 
9.55(-ll) 
2.87(-ll) 
3.91(-13) 
1.53(-12) 
1.02(-12) 
1.04(-16) 
1.19(-15) 
1.89(-9) 
6.76(-10) 
6.16(-10) 
1.88(-11) 
5.15(-11) 
4.42(-9) 
5.39(-ll) 



4.50(-l) 
4.26(-7) 
9.36(-8) 
1.34(-9) 
3.59(-4) 
1.68(-4) 
4.07(-5) 
3.29(-ll) 
1.87(-8) 
1.19(-10) 
1.15(-9) 
1.01(-12) 
3.41(-10) 
3.31(-12) 
6.39(-15) 
1.15(-14) 
5.51(-10) 
7.57(-10) 
8.10(-9) 
6.25(-10) 
3.09(-10) 
1.64(-9) 
7.43(-ll) 



1.96(-2) 
8.75(-8) 
8.06(-10) 
1.96(-9) 
9.91(-6) 
4.28(-6) 
7.27(-9) 
2.54(-10) 
4.97(-ll) 
6.64(-13) 
8.30(-16) 
2.03(-16) 
1.55(-17) 
1.80(-16) 
1.45(-23) 
5.74(-19) 
1.53(-11) 
4.77(-13) 
3.47(-13) 
3.09(-15) 
7.68(-17) 
3.77(-8) 
1.46(-13) 



4.99(-l) 
5.66(-9) 
1.64(-9) 
1.53(-12) 
3.14(-6) 
8.49(-5) 
5.40(-5) 
5.50(-13) 
1.40(-9) 
1.86(-14) 
2.69(-10) 
1.71(-13) 
5.20(-ll) 
2.60(-12) 
8.40(-14) 
2.09(-14) 
2.45(-ll) 
2.79(-ll) 
2.23(-9) 
2.42(-10) 
1.82(-10) 
9.93(-13) 
1.80(-14) 



4.99(-l) 
3.08(-9) 
1.44(-9) 
1.37(-12) 
4.03(-6) 
1.33(-4) 
4.32(-4) 
2.62(-13) 
1.01(-9) 
2.93(-14) 
2.09(-10) 
9.29(-14) 
1.13(-10) 
2.39(-ll) 
9.56(-14) 
8.42(-14) 
6.90(-12) 
2.02(-ll) 
1.82(-9) 
7.91(-10) 
5.54(-10) 
9.41(-13) 
2.49(-14) 



7.21(-14) 1.23(-12) 1.709(-18) 4.65(-12) 2.48(-12) 



4.94(-l) 
7.20(-8) 
2.55(-8) 
2.41(-12) 
8.18(-7) 
1.08(-5) 
2.59(-6) 
5.03(-13) 
3.30(-9) 
6.04(-14) 
7.28(-10) 
1.09(-10) 
2.16(-9) 
4.21(-12) 
3.77(-12) 
3.86(-14) 
1.27(-10) 
5.26(-ll) 
4.55(-9) 
1.06(-10) 
3.36(-10) 
1.08(-12) 
5.95(-14) 
1.43(-10) 



3.38(-l) 
3.63(-6) 
6.12(-7) 
7.36(-9) 
5.74(-5) 
1.62(-5) 
6.84(-5) 
5.98(-12) 
1.16(-9) 
1.83(-10) 
4.73(-ll) 
4.27(-12) 
7.01(-12) 
1.27(-11) 
5.64(-16) 
1.67(-15) 
6.79(-9) 
3.73(-9) 
1.20(-9) 
9.07(-ll) 
2.28(-10) 
5.34(-9) 
1.30(-10) 
2.08(-13) 



4.53(-l) 
6.61(-7) 
1.65(-7) 
1.78(-9) 
3.18(-4) 
1.24(-4) 
1.26(-4) 
3.14(-11) 
2.35(-8) 
1.75(-10) 
1.67(-9) 
2.78(-12) 
1.63(-9) 
7.34(-12) 
1.48(-14) 
1.53(-14) 
9.29(-10) 
1.83(-9) 
1.06(-8) 
1.59(-9) 
1.89(-10) 
1.92(-9) 
1.28(-10) 
1.87(-12) 



4.95(-l) 
1.37(-8) 
2.80(-9) 
2.93(-12) 
9.94(-7) 
1.27(-5) 
4.99(-7) 
5.13(-13) 
5.88(-10) 
1.18(-14) 
3.15(-11) 
7.18(-14) 
1.39(-12) 
1.85(-14) 
8.14(-15) 
1.34(-15) 
1.24(-10) 
1.05(-11) 
2.42(-10) 
2.07(-12) 
2.31(-12) 
1.78(-12) 
7.98(-15) 
9.33(-13) 



4.99(-l) 
7.53(-8) 
1.94(-6) 
1.85(-12) 
1.66(-8) 
6.85(-7) 
1.41(-4) 
8.94(-15) 
3.34(-10) 
4.97(-14) 
2.12(-12) 
7.93(-9) 
3.55(-7) 
4.70(-8) 
7.68(-ll) 
1.20(-11) 
8.97(-ll) 
1.64(-10) 
7.26(-9) 
2.57(-9) 
1.87(-9) 
1.12(-12) 
1.10(-13) 
1.22(-11) 



4.99(-l) 
3.82(-8) 
5.75(-6) 
2.69(-12) 
5.03(-9) 
9. 11 (-7) 
5.63(-4) 
4.67(-15) 
2.77(-ll) 
7.49(-14) 
8.35(-14) 
4.27(-9) 
1.37(-6) 
1.15(-7) 
3.06(-ll) 
1.86(-11) 
7.61(-11) 
1.16(-10) 
3.28(-9) 
8.64(-9) 
2.13(-9) 
1.68(-12) 
1.12(-13) 
1.32(-12) 
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Figure 2. Fractional abundances (n(X)/npj) of species X in logarithmic scale where njj is the total number of hydrogen atoms, obtained 
for the models with a metallicity of 0.1 zq for an A v of 3 mag as functions of cosmic ray ionisation rate (£) in s _1 (see text in Sect.[4]|. 
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Figure 3. Fractional abundances of species X obtained for the models wl" ' "' 

|4]and the caption of Fig. [JJ. 
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metallicity of 0.1 zq for an A v of 8 mag (see text in Sect. 
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Figure 4. Fractional abundances of species X obtained for the models with a metallicity of 0.1 zq for an A v of 20 mag (see text in Sect. 
[4] and the caption of Fig.[2]l. 



